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QUASI-STATIC HYDRODYNAMIC LIMITS 


ANNA DE MASI AND STEFANO OLLA 

Abstract. We consider hydrodynamic limits of interacting parti¬ 
cles systems with open boundaries, where the exterior parameters 
change in a time scale slower than the typical relaxation time scale. 
The limit deterministic profiles evolve quasi-static ally. These lim¬ 
its define rigorously the thermodynamic quasi static transforma¬ 
tions also for transition between non-equilibrium stationary states. 
We study first the case of the symmetric simple exclusion, where 
duality can be used, and then we use relative entropy methods to 
extend to other models like zero range systems. Finally we con¬ 
sider a chain of anharmonic oscillators in contact with a thermal 
Langevin bath with a temperature gradient and a slowly varying 
tension applied to one end. 


1. Introduction 

Quasi-static transformations are defined in thermodynamic litera¬ 
ture as those transformations where the external conditions change so 
slowly that at any moment the system is approximately in equilibrium. 
They are usually presented as idealization of real thermodynamic trans¬ 
formations [3j . On the other hand they are necessary, in the traditional 
approach, in order to construct thermodynamic potentials, for example 
to define thermodynamic entropy from Carnot cycles. 

Equilibrium statistical mechanics succesfully describes the equilib¬ 
rium thermodynamics states from the microscopic dynamical proper¬ 
ties of the system, in particular the corresponding stationary proba¬ 
bility distributions, given explicitely by the Gibbs measures. Let us 
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consider the simple situation when the internal dynamics of the sys¬ 
tem has a certain number of conserved quantities (energy, density,...), 
and the external conditions (temperature of the heat bath, pressure 
of the external forces,...) select the stationary equilibrium probability 
distribution for the microscopic dynamics. The fact that, when the 
microscopic system is large , there are only few of such conserved quan¬ 
tities, and correspondingly few external parameters will determine the 
equilibrium state, is related to the ergodicity of the microscopic dy¬ 
namics. 

The system can be driven out of equilibrium by changing the exterior 
parameters. If this is done abruptly, eventually the system will relax 
to a new equilibrium state corresponding to the new values of temper¬ 
ature, pressure etc.. The typical relaxation time to equilibrium may 
depend on the size of the microscopic system. 

The situation is similar but more complex when the external agents 
are space inhomogeneous, for example when the system is submitted to 
gradients of temperature or pressure. In the stationary situation fluxes 
of the conserved quantities can cross the system. The correspond¬ 
ing probability stationary distributions for the microscopic dynamics 
are called non-equilibrium stationary states. These are not explicitly 
computable as in the equilibrium case, but are locally close to the equi¬ 
librium ones, when the size of the system is very large. 

In the present article we want to obtain, at least for some systems, in 
a rigorous way the quasi-static thermodynamic transformations from 
the microscopic dynamics, through a space-time scale limit. In sys¬ 
tems of size N with diffusive behavior, where the typical time scale 
for reaching equilibrium is of order N 2 , we look at the evolution of 
the conserved quantities at the time scale N 2+a t, for a a > 0, and we 
change the exterior conditions smoothly at this time scale (i.e. very 
slowly with respect to the typical relaxation time). We obtain that in 
this time scale, at each instant t, the profile of the conserved quantities 
converges to the deterministic profiles corresponding to the stationary 
state for the boundary condition given at time t. This means that 
the profiles evolve quasi-statically, following the rules of the thermody¬ 
namic quasi-static transformations. Consequently the paradox of the 
quasi-static transformations disappears when these transformations are 
understood (also in a mathematically precise way) in a proper macro¬ 
scopic space-time scale. Remarkably we also obtain quasi-static trans¬ 
formations between non-equilibrium stationary states, that could be 
the base for constructing a thermodynamics of non-equilibrium. 

We consider three different models, in order to illustrate the gener¬ 
ality of the idea. 
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We consider first the symmetric simple exclusion in an interval of 
size N, i.e. particles performing symmetric random walks with ex¬ 
clusion rule, with density reservoirs at the boundaries. These reser¬ 
voirs are realized with random creation and annhiliation of particles so 
that the average density of particles at the boundaries will be respec¬ 
tively p-(t) on the left and p+(t) on the right. The rate of jumps or 
of creations/annhiliation is N 2+a , a being any positive number. We 
prove that at every time t the empirical particle’s density converges, as 
N —* oo, to the linear interpolation between p_(£) and p+(t), i.e. to the 
solution of the Laplace equation d xx p = 0 with the given boundary con¬ 
ditions. We actually prove the strongest result of the local equilibrium 
property. 

We first prove the above convergence in section [2] by using duality 
techniques as in [8], where analogous results were proved in the case of 
p± fixed constants. With the duality techniques we prove local equi¬ 
librium in the strong form of Theorem 12.11 and the same techniques 
allow to study the limit of the covariance function. In fact we prove 
in Theorem 12.61 that the covariance behaves at each time t as in the 
corresponding stationary states, as first computed by Spohn [15] . 

In Section [31 see Theorem 13.11 we prove a weaker version of the 
local equilibrium by using entropy methods as in [9] and m that are 
suitable to be used for more general models. In particular the proof 
of Theorem 13.11 can be adapted to prove the analogous statement for 
all gradient conservative dynamics with creations/annihilations at the 
boundaries as the ones studied in [7] in the case of constants p±. 

In Section [H we consider the zero-range process, where more particles 
per site are allowed, and jump symmetrically with rate depending on 
the number of particles in the site they occupy. A chemical potential A 
characterizes the equilibrium measures, which is the expectation of the 
jump rate, and in general is a non-linear function of the density. At the 
boundaries the reservoirs are given again by creations/annihilations of 
particles corresponding to two possibly different values of A ±(t) at time 
t. Jumps and creations/annihilations happen at rate proportional to 
N 2+a , and we prove that at each macroscopic time t (in a weak sense) 
the empirical density of the particles converges to the solution of the 
non-linear Laplace equation 

d xx X(p(x,t)) = 0, A(p(±l, t)) = A ±(t). (1.1) 

Finally, in Section [5l we consider a chain of (unpinned) anharmonic 
oscillators in contact with a heat bath with a gradient of temperature 
and subject to a force (tension) f (£) acting on the last right particle, 
that changes at the macroscopic time scale. The left end of the chain 
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is attached to a point and the heat bath is modeled by Langevin ther¬ 
mostats acting on each particles. The temperatures of the Langevin 
thermostats may depend on the atom and change slowly form one atom 
to the neighbor one, giving rise to a smooth macroscopic profile of tem¬ 
perature. When temperature profile / 5 _1 (r/) is constant in space and 
time at a value /3 _1 and tension is constant f(t) = r, the stationary 
(equilibrium) distribution is given by the canonical Gibbs measure with 
these values of temperaure and tension. In equilibrium the tension r 
is equal to the expectation value of the anharmonic force between par¬ 
ticles. 

For a non-constant profile of temperature we have a non-equilibrium 
stationary state where energy flows from hot to cold thermostats. We 
consider only situations where the profile of temperature does not 
change in time, and only the applied tension is changing in the macro¬ 
scopic time scale. Again the dynamics is speeded up in time by N 2+a . 
We prove that the empirical strain of volume converge to the solution 
of the equation 

r{r(x, t), P(x)) = f{t) (1.2) 

where r(r, /?) is the equilibrium tension corresponding to the volume 
(lenght) r and temperature /3 _1 , i.e. r = <9 r .F(r, j3), where T(r, (3) is 
the thermodynamic free energy. We have obtained this way the quasi¬ 
static isothermal transformation. 

The main interest in this last model is that we can define the heat 
Q(t) as the limit of the the total flux of energy between the system and 
the thermostats, divided by the number of particles, and it turns out 
that this is a deterministic function of time that satisfy 

U(t)-U(0) = Q(t) + W(t) (1.3) 

where U(t) is the limit of the internal energy (per particle) and W(£) is 
the work done by the exterior force. In the context of the equilibrium 
case (constant temperature profile), this is the first principle of thermo¬ 
dynamics for isothermal quasi-static transformations. For non-constant 
temperature profile Q(t) is usually called excess heat. Also this quasi¬ 
static transformation satisfies Clausius identity, i.e. the work W(f) is 
equal to the time difference of the free energy after properly identifying 
the macroscopic free energy as integral over space of the equilibrium 
one. 

Quasi-static limits where previously considered on the macroscopic 
hydrodynamic diffusive equation, rescaling the boundaries parameters 
after having performed the hydrodynamic limits (cf. [j], p4], [13]), 
i.e. in a two step approach. The point of this article is to show that 
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quasi-static transformations can be obtained in a straightforward way 
from the microscopic dynamics by considering the proper space-time 
scaling. 


2. Simple exclusion with boundaries 

We consider the exclusion process in {0, 1} Ajv , A ^ '■= {—N,..,N} 
with reservoirs at the boundaries with density p±(t) G (0,1). We as¬ 
sume that p±{t) are Lipschitz continuous. Presumably our results can 
be generalized to piecewise continuous functions but we did not inves¬ 
tigate this case. 

Denoting by rj(x) G {0,1} the occupation number at x G A-n we 
define the dynamics via the generator 

Ln,i — N 2+a [L exc + Lbj], t > 0, a > 0 (2-1) 

where 

- N—l N -1 

WM = 2 E (/(’U'+U - /O)) =: 2 E Ws-H m (2.2) 

x=—N x=—N 

r]( x ’ y } is the configuration obtained from r] by exchanging the occupation 
numbers at x and y, and 

L b ,,/W = E<WU*" V) ( 1 - p„(())’ ( " K )(/(," K ) - /(,)] (2.3) 

cr==b 

where rf(x) = 1 — rj(x), and r] x (y) = rj(y ) for x ^ y. 

Theorem 2.1. For any a > 0 and for any macroscopic time t > 0 the 
following holds. For any initial configuration rjo, for any r G [—1,1] 
and for any local function <p 

=< <p(v) >p(r,t)=: ^(p(d^)) (2.4) 

where [•] denotes integer part, 6 is the shift operator, < ■ > p is the 
expectation with respect to the product Bernoulli measure of density p, 
and 

PM) = \[p+(t) ~ P-{t)]r + ^ [p+{t) +p-(t )], re [-1,1] (2.5) 

is the quasi-static profile of density at time t. 

We prove Theorem 12. ll bv using the self-duality of the exclusion that 
can be adapted to our process as we explain in Propositions 12.21 and 
12.41 below. The dual process is defined in terms of absorbed stirring 
walks as explained in Definition 12.11 below. 






6 


ANNA DE MASI AND STEFANO OLLA 


In the rest of this section we use the equivalent definition of the 
dynamics obtained by considering the generator 


( 2 . 6 ) 



and then studying the process up to times of order N 2+a t } t > 0 being 
the macroscopic time. 

Definition 2.1. The absorbed walkers. 

One absorbed random walk is the Markov process in 
Atv+i := [-N — 1, N + 1] fl Z that after an exponential time of mean 
1 jumps with equal probability on its nearest neighbor sites and when 
it reaches ±(N + 1) it stays there forever. We denote by {x t ,t > 0} a 
trajectory and we denote by r the absorbing time: 

r = min{r_(jv+i), r N+1 }, r a = inf{f > 0 : x t = a} (2.7) 

For any positive integer n we also consider the process .,x^), 

t > 0} of n stirring walks absorbed at ±(iV + 1). The stirring process 
is defined by exchanging independently at rate \ the “content” in x and 

x + 1 , so that if xf* = x before the exchange then it moves to x + 1 
while if Xt — x + 1 then it moves to x after the exchange (see (I2.24p 
below for the generator of the stirring process). If one disregards labels 
this is the exclusion process. We denote by s > 0}, i = l,.,n the 
trajectory of particle i starting from x\) } = Xi and call Ti its absorption 
time. We call C the generator , Vx the law and Ex the expectation of 
this process starting from x = (aq, .,x n ). 

Proposition 2.2. Let r /o be any initial configuration. Then for all 
x G Ajv 


E vo [r]t(x)] = S x \rio(xt)l T> t] + £ x [l r <t u Xr (t - r)] ( 2 . 8 ) 


where 


'“±(«+i)M = (2+ " ) s) 


(2.9) 


Below we often write u±(s) = m±(n+i)(s)- 

Proof. Given x and t, a trajectory x = {x s ,s 6 [0, t]} of the 
absorbed random walk starting at x and a trajectory 77 = {r] s , s E [0, t]} 
starting from 770 , we define for all s < t 



( 2 . 10 ) 
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Calling P = V x x P % the (product) law of (x, 17 ) and E the expectation 
with respect to P, we observe that for x G Ajv 


= E r;o (rj t (x)), E(iJj(x,r), 0)) = right hand side of 

( 2 . 11 ) 

Thus to prove (12.81) we show below that j^E [(^(x, 17 , s))] = 0. We 
first define for x G Ajv+i, r] G {0, 1} Ajv and a± G (0,1) 


<Pa±(x,v) = 


\v(x), 

VI 

^4-H 

\a± 

if x = ±(N + 1) 


and we will prove that 

d_ 

ds 

Let 0 < s < s' < t and call 


E{ip{K,V,s)) = ( ^ E(ip a± (x t - S ,ris )) 


a±=u±(t—s) 


( 2 . 12 ) 


(2.13) 


D(s.s') = ip(x,rj,s') - ili(x,rj,s), 

D(s,s') = - <Pa ± (x t (2.14) 

First observe that for any choice of a±, 

D(s, s') W s ' = 0, D(s, s')l T <t—s' = 0 (2.15) 

If instead t > t — s' it is important to choose a± — u±(t — s). In fact, 
since u± are Lipschitz continuous functions, we have 

sup ^(x^cr) ~ < -Pu±{t-s)(xt-<nr) (J )\l T> t- s ' < C{s' - s)l Te [ t _ s /, t _ s] 

o-G[s,s') 

(2.16) 

Since P(t G [t — s',t — s]) — > 0 as s' 4- s, from from (I2.15p and (12.161) 
we get that 


, lm U£M) 

s'ps s' — s 



E(D(s,s')) \ 
s' - s J 


a±=u±(t—s) 


Thus the Proposition will follows from the fact that the derivative of 
E(ip a± ) is 0. This last derivative is easier to compute because Lp a± is a 
function of (x, rj) unlike 7/1 that is a function of the whole trajectories 
(x, rf). In fact we have 

-^-E[if a± (x t -s,Vs)\ = -E[Cip a± (x t s,r]s)\ + E[L NtS cp a± (x t - s ,r] s )] 

= -[a+ - E(r) s (N)l Xt _ B=N )} + [u + (t - s) - E(r] s (N)l Xt _ s=N )] 
~[a~ - E(rj s (-N)l Xt _ s= _ N )] + [u_(t - s) - E(r] s (-N)l Xt _ s= _ N )} 

(2.17) 
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because the contribution of the generator L exc cancels with the jumps 
inside [— N, N] of the random walk. Thus (j^E[<p a± (x t _ s , r} s )] = 0 by 
putting in (12.171) a± = u±{t — s ) in agreement with (12.131) . □ 


Corollary 2.3. For any a > 0, any initial configuration rjo, any t > 0 
and any r G [—1,1] 

lim E V0 [r] N 2 + at ([Nr})] = p(r,t ) (2.18) 

N—t oo 

Proof. We use (12.8[) . First observe that for any x G 

V x {r >s)< ce~ c ' s/N2 (2.19) 

thus the first term on the right hand side of (12.81) converges to 0 expo¬ 
nentially. Call 


F n ,x(s) — V x [r N+1 < s, t_tv-i > r N+1 ] 


then by (12.191) 


\£x [lr=r JV+1 <iV2+-t u N+1 (N 2+a t - r)] 


rN 2+ %t 


' 0 


u N+l (N 2+a t - s) dF N , x (s)\ 


-c'N a / 2 


< ce 

Recalling (12.9[) we next observe that by Lipschitz continuity 

r N2+a/2t tN 2+a / 2 t 

/ u N+1 {N 2+ °t - s) - p + (t ) dF NyX (s) < c ——— < c 


_/y2+o - “_/Va/2 

Finally, by (12T9D . 

\F N , x {N 2+a / 2 t ) - V x [t_ ( jv+1) > TAT+r] | < V x [tjv+i > N 2+a ' 2 t] < ce~ c ' Na/2 
and analogously for the term with w_pv+i). Finally observe that 

X 

V x {t N+ i < T_(JV+1)) = - + 


2 2(N + 1) 


□ 


We now generalize (12. 8 j) by writing a duality formula for the corre¬ 
lation functions. We thus consider n > 1 distinct points x±,..,x n in A N 
and recalling Definition 12.11 we prove the following Proposition. 

Proposition 2.4. For any positive integer n, any x\,..,x n distinct 
points in An and for any go we have that for all t > 0 (recall (12.9j ) ) 

n n 

= 4[n + 1 Ti<t U x 0)(t - Ti )] 

i= 1 i=l 


( 2 . 20 ) 
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Proof. Denoting by E — Vx x E r;o and by X_ = (x^, x^\ s G [0, t]), 
we prove below that 


^EftiX^s)] = 0 , 


Z— 1 


( 2 . 21 ) 


The strategy is to adapt the proof given in Proposition 12.21 for n = 1. 
Let J r /_ s be the sigma-algebra generated by {(x^,, ,x^),s' e [0,t — 
s]}. Given we know k— number of particles absorbed by time t — s 
and the times of absorption r, < t — s, i = i±,..ij.■ Thus calling 


U(t - s) := \\_ u x U){t ~ Tj)l Tj <t- s , I = {k, -h} 
jei Tj 


we have that for all 0 < s' < s < t, 

E^X^s^-ViX^s))] 


= E 


U(t - s)E[Y[ 1 s') - s)) | Ji_ s ] 

i£I 


Recalling (12.1211 . as in (12. 16 j) we have that for all 0 < s' < s < t 
sup JJl Ti>t -s^{x {l \ria) - Y[l Ti >t-s<P u± (t-s)(xi-a, r l,c r ) 

o-e[ s ', s ) . . 

— C(s s) ~^3i:TiG[t—s,t—s'] 

Thus as in the case n = 1, (I2.20p is a consequence of 


lim — 

s'/'s s — s' 


-E 


U(t — s)E[ JJl T . >t _ s 

i<£I 

<-(i) 


(<Pu±(t-s)(x {l) (t - s'),v) - (pu±(t-s)(x {l) (t - s),77 8 ))|7i_J 


= 0 


The above equality can be proved by applying the generators similarly 
to (12T711 . ^ " □ 


Corollary 2.5. For any initial configuration r/o, for any t > 0, any 
n > 1 and any distinct points x\,..,x n in An, calling Lv = N 2+a t 


lim 

iV—>■ oo 


E, 0 [ Vt N (^i)] E, )0 


[vt N (xi)] 


i —1 


i=l 


= o 


( 2 . 22 ) 


Proof. By (12.2011 and (12.19j) it is enough to prove that 

lim I D x n I = 0 

TV—>• oo — ’ 


(2.23) 
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where 

n n 

Dx,N I 711 ^Xi ^" 2 )] 

2=1 2=1 

We use the coupling explained in [ 6 ] (see also 0 ) between n inde¬ 
pendent random walks and n stirring particles. Both the stirring and 
the independent particles move in the whole Z. The evolution of the n 
stirring particles is the process with generator L st that acts on functions 
/ on = {x E Z n : X{ ^ Xj for all i ^ j} as follows 

1 n 

L*m = 5 E E [/($*’“■) - /($)] (2-24) 

2=1 cr==bl 

with x l,er = x + aej — ere,,- if there is j ^ i such that Xj = Xi + ae^ 
otherwise x l,cr = x + <re;, e* being the unit vector in the i direction. 
We refer to Section 6.6 of [ 6 | for the dehnition of the coupling and we 
denote by x®{t) and Xi(t) the position at time t of the independent and 
respectively interacting particle with label i. The initial positions are 
the same x 0 ( 0 ) = x( 0 ) and the law of the coupling is denoted by Qx 
(we denote by Q% also the expectation). In Proposition 6.6.3 of [ 6 ] it 
is proved that for any positive 7 and k there is c so that for all t > 0 

Qx[ sup |xj(s) — x?(s)| < t* +ry , Vi = l,..,n] > 1 — ct~ k (2.25) 
\0<s<t J 

We call Tj 0 = min(r°’ _ , r l 0,+ ) with r ? 0,± the first time the independent 
particle i is at ±(N + 1), r* = min(r i + ,Tj _ ) is the analogue hitting 
time for the interacting particle i. Observe that since the particles are 
moving by stirring, then the hitting time of interacting particle i is not 
influenced by the fact that another particle (say j ) is inside the interval 
after being outside. Moreover the function inside the expectation in 
Dx,n depends only on the hitting times. Thus 

n 

Dx,N — Qx(^ Y J_ U Xi(ri){tN ~ T~i )] 

2=1 

n 

- II U x°(T?)(t N - T?)]) 

2=1 

We call a — (a\,a n ), at E {—1,1} the vector that specihes where the 
stirring particles exit the interval, namely a* = ±1 if stirring particle 
i exits in ±(N + 1), thus t* = r“‘ . Analogously b = (bi,..,b n ), bi E 
{—1,1} specihes that r 4 ° = r, 0, ’ . Given a and b we denote by r(a) = 
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( t “\ r“ n ) and by r( 6 ) = (rf 1 , -,r bn ). With this notation we have 

= E E <3*(0(t(o), r ©)) ( 2 - 2 6 ) 


a 6 


i=1 


&(T(a),T(b)) = [1 Ti=T °i< tN U ai (t N - Ti)\ ~ 

n 

II 1 1 r?=r°<t jv Ub ^ N ~ T i 


i=l 


Call 


A n = {n < X 2+ "/ 2 , r° < N 2+a ' 2 for all i = l,..,n} 

Then by (12.191) 

Qx(<i>(T(a), r(b))) = Qx.(^(r(a), r(b))l AN ) + c'e _cJV “ /2 (2.27) 
We now use the Lipschitz continuity of u± 

\Qx{(/>( r (a),T(b))l AN ) - Qx{$(r(a),T(b))l AN ) \ < c"^-^ 

n n 

®(r(a),T(b)) = Y[u ai (t N )l T . =T H< tN - IJw 6 i (Cv)l r o =r o, 


i— 1 


i —1 


■x=x i <^ 


Observe now that if a — b then < h(r(a), r(b)) = 0, thus 


_ —b _ _ 1 0_ 0,b_ _ 

= T i rP= r i <* 


\Qx(${r(a),T(b))l AN ) < OE E Qmi 1 * 

i= 1 6e{—1,1} 

Call 

= {|xj(f) — x°(f)| < X^+f for all f E [0, iV 2+CK//2 ]} 
thus by (I2.25p for a suitable 7 and for /c sufficiently large 


J 

(2.28) 


E ^( 1 r,;=r- 6 <Lv 1 T0=r°’ b <t JV ) 


66{-l,l} 


-Qx{1b n E i 1 
&€{-!, 1 } 


T»=n b <ijv 


T?=T?' b <t, 


, 1 ) 


< cX“ fc 
(2.29) 


Call X± = TV + 1 t A^ 2 + f and Y± = -X - 1 t Xs+f. Call 
CV = { independent particle i reaches and goes to X^ before reaching Yjj } 
{J {independent particle i reaches X^ and goes to before reaching A"^ j . 
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Then in C%, the complement set of Cm, we have that independent and 
interacting particles i exit from the same side. Thus 

Qx{f-C^B N ^ ( ^-Ti=rP<t N ) — Qx{1cc n 1b n y ] ^T i =rj'<tj V lrP=r I 0 ’ i '<tjv') 
— 1 , 1 } &£{— 1 , 1 } 

so that 

Qx{1c%1b n [ 1 r i=r~ b <t N ~ 1 rO=r°’ b <l Ar ]) = 0 

6 e{-i,i} 


We are thus left with the trajectories in Cm which is an event that 
involves only the independent particle i. Thus 

Qx(1 Cn 1 Bn [K i= rr b <t N -Q^ Cn ) 

be{~ 1,1} 




jVs+f 
2(N + 1) 


□ 


Proof of Theorem ld.il Since any cylinder function is a linear combi¬ 
nation of functions of the type r)(x), (12.4j) is a consequence of (I2.18j) 
and (12.22)1 . □ 


We conclude this section, by computing at any macroscopic time the 
covariance of the process. For xi ^ x 2 E A N and s > 0 we call 

2 

v n (x 1 ,x 2 ;s) :=E. w ^}J [r/ s (x i )-p N (x i ,s)]'), p N (xi,s) — E r/0 (rj s (x i )) 

i=1 

(2.30) 

We know from Corollary 12.51 that vm converges to 0 as N —» oo and in 
the next Theorem we show that it is of order IV -1 . 

Theorem 2.6. For any initial configuration r/o, any n < r 2 E (—1,1) 
and any t > 0 we have 

lim^NvM ([Ain], [Nr 2 ]; N 2+a t) = ~ [p + (i)-p_(i)] 2 ri(l-r 2 ) (2.31) 
Furthermore for any t > s 

lim NK V0 ([r] N 2 +as (x)-p N (x,N 2+a s)}[r] N 2 + at (y)-p N (y,N 2+a t)\) =0 

(2.32) 
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Proof. We compute the /.-derivative of vn- 

= ^vo{L N , t [Y[(v(xi,t) -p(zi,t))]) 

i=1 

2 

- ^ Er /0 [Vfa, t) ~ p(Xj, t)] 1 j^i) 

*j'=l 


After easy calculations we get: 


dv N (x\, x 2 ; t) 
dt 


L st:N v N (x 1 ,x 2 ]t) - ~[p N (xi,t ) - Pv(at 2 ,t)] * 2 l|*i- 3:2 |=i 
-ujv(ati, x 2 ; t) [l|xi|=iv + l|x 2 |=iv] (2.33) 


where vn{xi,x 2 ] t) is thought of as a function of the positions or and x 2 

and L st,N is the generator defined in (12.24)1 but with “jumps” outside 

2 

Ajy suppressed. By setting vn(xi, x 2l t) = 0 in the set UlN = JV + i} 

i=l 

we can rewrite (12.33P as an equation for vn(xi, x 2 ,t) for x\ ^ x 2 with 
\xi\ < N + 1 as follows 


d 

dt 


v N (xi,x 2 ]t) = 


v N {xi ,x 2 ,t) = 0 


I 2 

L st v N (xi,x 2 ,t) - -[p N (x u t) - piv(a: 2 ,t)]''l|* 1 -* 2 |=i 
2 

in Uiu i| = N + 1} and when t — 0 (2.34) 

i=l 


which is solved by 


v N (xi,x 2 ,t) 


E x llX2 ( [pn(xi{s), t - s) - p N (x 2 (s),t - s)]' 


xl 


|xi(s) — a;2 (s) 1=1 J-ri>S,T 2 >S 


(2.35) 


where Ef y denotes expectation with respect to the process of two stir¬ 
ring particles, labelled 1 and 2, which start from x\ ^ x 2 E An respec¬ 
tively. T\ and r 2 are the first time when particle 1, respectively particle 
2, reaches ±(N + 1). From (12.181) we get that for = N 2+a t 


lim sup 

N ^°° s <N 2 + a / 2 t 


sup N 2 

xgAjv 


[p N (x,t N - s) - p N (x +l,t N - s)]‘ 


1 

2 


[p+{t) ~P-(t)] 



(2.36) 
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From (12.191) we get 

•t N 

E S 1 [p7v(xi(s),t - s) - p N (x 2 (s),t- s)]' i l|x 1 ( a )-a 2 (a)|=ll 


XI ,X2 ' 

fl jY 2 +a /2 


n>s,72>s 


nst 


'0 


[pN^X^s)^ - s) - p N (x 2 (s),t- S)]' 


^ l|xi(s)—X2(s)|=lln>s,T2>s 


, / — cN a / 2 

< ce 


(2.37) 


From (12.361) and (I2.37I) . 

lim AT;at ([ATy], \Nr 2 ];N 2+a t) 

N^roo v 7 

^ roo . . 

= — ^ [p+(t) — P-(t)] J E rij7 - 2 y§B(s) {B (s)) 1 t >s,t'>s J 

Sb(x) the Dirac delta at x, B(s ) and B'(s) independent Brownian mo¬ 
tions starting from 74 and r 2 and r, r' the hitting times at ±1. As 
the right hand side of the above expression is identify to the kernel of 
(—A) -1 with Dirichlet boundary conditions, (12.31(1 easily follows. 

To prove (12.32(1 we observe that 


E r,o ( [vn2+« s (x) - Pn(x, 1V 2+ "s)] [r] N 2 + at (y) - p N (y , N 2+a t)] ) 

= E vo (\r] N 2 + as (x) - p N (x, N 2+a s)]p N 2 + at (y)^j 

= E vo ( [rj N 2 + as (x ) - p N (x, N 2+a s)]£ y (l T>N 2 + a{t _ s) r] N 2 + as (x(N 2+a (t - s)) 
TTy (lr<A 2 + a (t— s)'U j Xt if f) )Ejj 0 ^ [pA r2 +“s (x) PN^Xj s)]^ 

= E w (jrj N 2 +a s (x) - p N (x, N 2+a s)}£ y (l T>N 2 + a{t _ s) p N 2 + as (x(N 2+a (t - s)) 


that by (12.191) implies (I2.32|) . □ 

As a corollary of Theorem 12.61 one can prove that the fluctuation field 
converges to a Gaussian field as in the case of constants p± studied in 
[12] • We do not proof this here. 


3. Simple exclusion: entropy method 

We prove here the quasi-static limit using the entropy method. This 
gives a weaker result, but the methods extend to other models where 
duality cannot be used. We assume here that p±{t ) are differentiable 
functions of time with bounded derivative. This can be relaxed to 
Lipschitz continuity, as assumed in the previous section, but the proof 
would be more involved. 

In the following rj t is the process generated by (12. ip . 
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Theorem 3.1. For any a > 0, any t > 0, and any local function (p(r]), 


lim E. 

TV—>• oo 


'VO 


ds 


N 




G{y)<P{p{y,s))dy 


where G is a continuous test function on [—1,1]. 


= 0 . 
(3.1) 


The above theorem is stated for any initial configuration r/o. More 
generally we can start with any initial distribution, since in the follow¬ 
ing all we need is that the relative entropies are bounded by CN for 
some constant N, that in this case is automatically satisfied by any 
probability measure on the configuration space {0, l} 2Ar+1 . 

Consider the empirical distribution of the density at time t: 

1 N 

^ (g > = 2VTT &')*<*> (3-2) 

x=—N 

where G(y, t ) is a smooth function on [—1, l]xR + with compact support 
in (—1,1). The time evolution is given by 

Ut(G)-Uo(G)= f^N^+d^ds+OiN^+M^t) (3.3) 

Jo 

where Mjv(t) is a bounded martingale. By dividing by N a , we have 
immediately that 

lim [ £ Ns (dyG) ds = 0. (3.4) 

N ^°° Jo 

This implies that every limit point of the empirical distribution £jv,t is 
a measure on [—1,1] that satisfies Laplace equation, in a weak sense, 
with boundary conditions that we will identify in the following. 

Corresponding to the quasi-stationary profile p(x,t ) defined in (12.5ft . 
we consider the inhomogeneous product measure 

N 

M»)= n <X-LT ( U1 -p(|f,«)] 1 -’ ( * ) (3.5) 

x=—N 

as reference measure at time t. The Dirichlet forms associated to the 
generator are then l D± N ,t,p±(t)(f) and T> e;E ,t(/): 

».,<,„(/) = i^y-’>W( 1 _p)’>W[ v / 7M - y/Mf/U-)), x = ±N 

V 

1 Ar_1 2 
®ex,t(f) = ( v ^+iV / 7M) fhiv) 

Tj x= — N 


(3.6) 
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Let /jv,t be such that the law of rj t is given by /tv,* Ht, the following 
holds. 


Proposition 3.1. There is C so that for all t 

I N,s,p+(s) (In, s') “1 ® — N,s,p— (s) (/)v,s) d" ^ ex,s(^fN, s)^) ^ 


ct 

iV' 


(3.7) 


Proof. We have that 

d t ( f N ,t lh) = (■ L *N,tf N ,t ) Pt (3.8) 

where is the adjoint with respect to /x f (r/) that can be computed 


as: 


£<?(>*)( W)( r,Mv) = Y. F (v)( L N.,G)(riMri) 

V V 

+JV 2 +” £ Ffo) £ G(ri x ’ x+1 ) ' I+1 ) - /*(»>)) 

77 tc 

Observe that 






)■ 


- 1 


Mv) 


and 

P(“iv~" j ^)(1 ~~ P(jv> 0 \v{ x )— 17 ( 3 + 1 ) 

P(#» 0 (! — 0 


) ? 


1=[1 + 


p'(jM) 


iVp(f,f)(l-p(§,f) 

p'(IM) 


)] 


? 7 (x)—??(a;+l) 


= ^W-^ + i)]. ( , t)(1 J- (f!i))+ o 


1 

iW 




Denote Hjv(t) = E ?/ /iv.tfa) log/ a^(??)Mp) 

j t H N(t) = log f N ,t(v)(LN, t f N ,t)(v)»t(v) + 53 [9t/jv,t(^)]^(7?) 

V V 

= ^f N A r i)( LN ,ti°gf N ,t)(v)vt(v) + t(p)K(p) 


Observe that, since ^ fN,t(v)pt(v) — 0; we Lave 

^[di./W,t(7/)K(p) = = 0(N). 
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By the inequality alog(6/a) < 2y/a(Vb — y/a), we have that 

f (v) (L N ,t log /) (rj) < 2v / 7M(^iv,iV / 7)(^)- 
Furthermore 

Y 2 \/ faAv) ( L N ,t y/f N ,t) (v)vt (v) 

v 

N-1 

— _ N 2+a E E V x } x-|-i 

rj x=—N 

- N 2+a (SjV, i]P+ (f)(/(V,() + l £>-N,t,p^(t){fN,t)) 

and 

Af-l JV-1 

EE EE (v a;,a:+l \/ /lV,t j Ah 

77 x=—N 77 x=—N 

N—l 

+E E a/ /Ar,f(h' C ’ X+1 )V a;>3;+ i(v / fN,t)(vNx,x+l^t(v)) 

rj x=—N 

N-l 

= EE (v #,#+1 a// lV,t) hi 

»? x=—N 

N-l 

+vE E V7U(T’ I+1 )v« + i(v / 7U)hW-^ + i)]s(^.«U< + o(iv- 1 ) 

rj x= — N 

= 2 ®ex,t(fN,t) + ^B N (t) + OiN- 1 ) 

with 

iv -1 

B N (t) = Y Y Vf^t(v x ’ x+1 )^x,x+i(^jN^t)(v(x)-v(x + l))B(^,t)n t 

rj x=—N 

By an elementary inequality we have: 


JV-l 


5iv(t) < ^®exAfN,t) + 7r^YYl fMV X,X+1 ) (v(x) - v(x + !)) 2 5 (t7, t) 2 Vt 


rj x=—N 


N-l 




x=—N 


N-l 


<(/«) + ^ E 1 i + E/«M v ™+ift 


x=-N 
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and since B(x,t ) < ^||p'||oo? iterating on the bound for V XjX +i//t, we 
obtain 


B N [t) 


N 

< —®ex,t(/ N,t) + Ca 
CL 


All together we have 


H N {t) — Hn(0) — —N^ +a / ((SjV,s,p + (s)(/lV,s) + ®-JV,s,p_(s)(/jV,s) + tD ex,s(fN,s)) ds 

Jo 

-N 1+a [ B n (s) ds + 0(N 1+a )t 

Jo 

< —N 2+a f (2)jv,p + (s)(/at,s) + ®-jV,s,p_(s)(/?V,s) + 2} e x,s(/v,s)) ds 

Jo 

]u2+a ft 

+ -z- / Dex, S (/A, S )ds + iV 1+Q 2Ct + 0(iV 1+ “)t 



JV,p_|_(s) (./jV,s) T ^ —N,p-(s)(f N,s) T 2® err (/jv, s )J ds 

^ ^jv(O) + At 1+Q 2Ct + 0(iV 1+a )t ^ C" C"'t 

— j\r 2+ “ — jV 1+Q! At 


□ 


Proof of Theorem \3.1[ Consider for any point y G [—1 + k/N, 1 — k/N\ 
the set A y of 2k + 1 integers defined by A y = {[ Ny ] — k, ..., [At?/] + /c}, 
and define /?v,*| a„ the marginal of f^.t on A y . It follows from (13.7ft that 
for ft... = lim f m 11 a 


A:—1 

*,«.»(/«,») - E E 

r/e{0,l} 2fc+1 x=—fc 

(where p k p( ^ y ^ is the Bernoulli measure on {0, l } 2fc+1 with density p(y, t)). 
This implies that f tt y(r}-k, ■ ■ ■ ,3k) is symmetric for exchanges. Fur¬ 
thermore, considering the boundary blocks {—At,..., —N + 2k} and 
{At — 2k ,..., At} we obtain respectively 

25jV,i,p + (t) (/?,—l) = 0, 2)— N,t,p-(t) (/i,l) = 0 

that implies ffi and f t ~i are constant, and since they are probability 
densities with respect to p k p+ ^ and they are equal to 1 in those 

boundary blocks. 
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From the above argument we obtain that 

nt /»1— k/N 

iim lim ds d yY] fN,s(v)Ps(v) 

k^ooN^ooJ 0 J_ 1+k/N ^ 

9 Mv s ) ~ <P ( ^ ^ vfr) 

\x—Ny\<k y \x—Ny\<k 

Next we extend this statement to macroscopic blocks. Consider now, 
for any y, y' E [—1 + k/N, 1 — k/N] the two blocks A y and A y >, and 
let f N,t,y,y' {j]i fy) the corresponding joint marginal. Dehne, for functions 
on two separate blocks f(rj-k, ■ ■ ■ ,rjk] • • •, fjk), the Dirichlet form 
corresponding to the exchange of the occupation of the centers of the 
box: 

(/) = X] (\/7( T o(h^)) - Vf{v,v)) A k p(t, y )(v)A k p(t,y')(v) 

v,v 



where Tq is the exchange of rj( 0) with 7/(0). 

By the same telescoping argument used in HP , we have that ft, y , y ' (v, v) 
lim f Nt t, y , y '(Vi v)i satisfies 

®o(ft,y,y') < c\y - y' I 2 


It follows that (see [TO] for standard details): 

-1 -k/N 


lim lim lim 

[ ds 

e—>-0 k —^oo TV—>-oo t 

Jo 

/ i 

[Ny]+k 

1 

E 

1 2k + 1 

V 

= [Ny]~ 


'-1+k/N 


dy^fN, s (v)As{v) 


[JV(y+e)] 


(3.10) 


2Ne 


d( x ) 


= 0 


x=[N(y-e)] 


In fact we already know that 

t / 1 ~N+2k \ 2 

k'S.lToo/, *EhrrT E "M-P-M = 0 

77 \ x=-N / 

(3.11) 

and similarly for the block [At — 2k, At], By (13.lOj) we have 


lim lim 

e—>-0 iV—>-oo 




-N{ l+2e) 

v( x )~P-( s ) 

x=—N 


2 


(3.12) 
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This implies, with ( 13.4j) . that the limit of 7 ^ Y^}[N(y-e)] r l( x ) as N y 
00 and e —> 0 converges strongly to the solution of the Laplace equation 
with this boundary condition given by p±(t), i.e. p(y,t). 

□ 


4. Zero range 

The results of the previous section extend quite straightforwardly to 
any gradient conservative dynamic. As an example, let us consider the 
zero range process, whose generator is given by 

fj2+a , N 

L°N, t f(rj ) = ——{ 9 (v(z)) [f( 9 x,x+(r ) ~ f{v)\ 

X=—N < 7=±1 ^ ^ ^ 

+ £M0[/(<r + )-/wi} 

cr=d= 

where g x,y is the configuration 77 with a particle moved from x to y, 
for x, y = For r ) N,N+1 and r ]~ N ~ N ~ 1 we have destroyed 

a particle in the correspondig site, while rj x,+ means we have addeed 
a particle in the site x = N, —N. We assume that the rate function 
g : N — y M + such that p(0) = 0, g{k) > 0 for k > 0, and sup fc | g(k + 
1) — g{k) | < +00. As before we assume that A±(i) are differentiable 
functions of time with bounded derivative. 

For any A > 0 consider the measure /i\ on the non negative integer 

9 ^' = ^°) ! = L 

where Z( A) is the normalization constant. 

Let A (r,t), r e [— 1 , 1 ] be the linear interpolation of A±(i) defined as 

in m- 

Consider as reference measure the inhomogeneous product 

N 

v?(v) = n VMfrMx)) (4.2) 

x=—N 

Observe that for all N, t and all local function ip: 

= 0. (4.3) 

v 

For constant A, we denote by p\ the corresponding homogeneoous 
product measure. 

Assume that the initial configuration 770 is randomly distributed by 
a probability measure with density f N}0 with respect to (77) and such 
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that the relative entropy ^ fN,o(v) l°g fN,o(v)ho(v) — CN for some 
constant C > 0. 


Theorem 4.1. For any local function ip, denoting ip(\) =< <p >\, the 
average of ip with respect to p\: 


lim E 

TV—>-oo 


ds 


G (jf) d MVs) - J^ G(y)<p(\(y, s))dy 


= 0 . 


(4.4) 

where G is a measurable bounded test function and $i is the space shift 
by i, well defined for N large enough. 


We omit the proof of Theorem 14.11 being very similar to the one of 
Theorem 13.11 for simple exclusion. 

As observed in [0] in the case of A± constants, (14. 3 p straightforwardly 
yields (14.4p . 


5. Damped anharmonic chain in temperature gradient 

We consider a chain of N coupled oscillators in one dimension. Each 
particle has the same mass, equal to one. The configuration in the 
phase space is described by rj = {q x ,p x , x — 1,, N} e K 2A . The in¬ 
teraction between two particles x and a; — 1 is described by the potential 
energy V(q x — q x -i) of an anharmonic spring. The chain is attached 
on the left to a fixed point, so we set qo(t) = 0,po(t) = 0. We call 
{r x = q x ~ q x - i,x = 1,..., N} the interparticle distance. We assume 
V to be a positive smooth function, and that there exists a constant 
C > 0 such that: 

lim —j—j— = oo, lim V"{r) < C < oo (5-1) 

|r|—>oo H |r|—s>oo 

Energy is defined by the following Hamiltonian: 

N 

X=1 

The particle dynamics is subject to an interaction with an environ¬ 
ment given by Langevin heat bath at different temperatures ft- 1 . We 
choose /3 X as slowly varying on a macroscopic scale, i.e. j3 x = f3(x/N ) 
for a given smooth strictly positive function (3(y), y G [0,1]. 


Px 


+ V(r x ) 


(5.2) 
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The equations of motion are given by 


dr x (t) = N 2 +a (p x (t ) -p x _i(t))dt 

dp x (t ) = N 2 + a (V'(r x+ 1 (t)) - V'(r x (t))) dt 

-N 2 +a -yp x {t)dt + N 1+a / 2 \J^.dw x (t), x = 1, N - 1 
dp Nit) = N 2 +a (f(t) — y'(riv(i))) dt — N 2 +a ^p N it) dt 



(5.3) 

Here {w x (t)}x are ^-independent Wiener processes, 7 > 0 is the cou¬ 
pling parameter with the Langevin thermostats. The force f (t) is as¬ 
sumed to be a smooth function of the macroscopic time t. 

The generator of the process is given by 


£ f«) , = Jv2+a ^0 + y S N ) , (5.4) 


where A T N is the Liouville generator 


N N—l 

A t n = Y (Px -Px-i)dr N + ^(H'(r x+ i) - V'(r x ))d Px + (r - 

rr=l x=l 

(5.5) 

while is the operator 


N 

Sn = Y (Px' 9 p* - Px d P,) ( 5 - 6 ) 

X=1 


A11 equivalent dynamics is given by a different modelling of the heat 
bath: particle x undergoes stochastic elastic collisions with particles of 
the environment at temperature fd~ l , i.e. at independent exponentially 
distributed times of intensity 7 N 2+a particle x changes its velocity to 
a new random velocity normally distributed with variance /3” 1 . The 
evolution equations are given by: 


dr x {t ) = 

N 2+a (p x {t ) -p x _i(t))df 


dpxit) = 

N^(V'(r x+1 (t))-V'(r x {t)))dt 



+ {Px^NZ+m) - Px(t~)) dAf x { , yN 2+a t), 

x — 1,.., N — 1 

dp N (t) = 

N 2+a {f{t ) - V'{r N ))dt 
+ {Pn,Mn(' yN 2 +°‘t) -PN(t~)) dJ\[ N {^N 2+a t) 

(5.7) 


where p X) k are independent gaussian variables on mean zero and vari¬ 
ance /3' 1 , and {J\f x (t), x = 1,..., N} are independent Poisson processes 
of intensity 1. 






QUASI-STATIC LIMITS 


23 


For f{t) = r constant, and (3 X = f3 homogeneous, the system has a 
unique invariant measure given by a product of invariant Gibbs mea¬ 
sures ij* 8 : 

N 

dn^p = n e~^- Tr ^- g ^dr x dp x (5.8) 

X=1 

where £ x is the energy of the particle x: 

2 

£ x = j + v (r x ). (5.9) 

The function Q Tt p is the Gibbs potential defined as: 


Qr,p = log 


\J2,71 / 3 _ 1 / e-WV-^dr 


(5.10) 


The free energy of the equilibrium state (r, /3) is given by the Le¬ 
gendre transform of —fd~ 1 Q r ^: 

T rfi = sup {nr + /3 _i £ Ti(9 } (5.11) 

T 


The corresponding convex conjugate variables are the lenght 

r(r, /3) = 

and the tension 

r(r,/3) = d r T r fi. 

Observe that 

E <f [V l (r)\ = T. 


(5.12) 

(5.13) 

(5.14) 


Let Vp the inhomogeneous Gibbs measure 


dvp. 


N 

n 

x=l 


-> — 0X&X 


(5.15) 


Observe that this is not the stationary measure for this dynamics for 
f = 0. 

Also we will use the thermodynamic entropy defined as 


S(u,r) 


inf {- 

/3>0 L 


-Pu -£.F(r,/3)}. 


(5.16) 


Let f t N the density, with respect to Up , of the probability distribution 
of the system at time t, i.e. the solution of 



rN 

t-N Jt ) 


(5.17) 
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(5.18) 


where jCffi’* is the adjoint of with respect to , i.e. explicitely 

N—l 

{cT)* = -N 2+a A T f - N 2+a ^(/3 x+ i - p x ) Px V'(r x+1 ) 

X=1 

+N 2+a f3(l)p x f(t) + N 2+a ^S N . 

Define the relative entropy of f t N du^ with respect to dvp : 

H, v(t) = I f, K log f?d$. 

We assume that the initial density satisfy the bound 

H n ( 0) < CN. (5.20) 

We also need some regularity of dehne the hypercoercive Fisher 
information functional: 

A—1 „ / o p AT . o pl\T\ 2 


(5.19) 


Mt) = Y,p ; 1 f 

X=1 ' 


(Aj?+d,jry 


f, 


N 


■dvp. 


(5.21) 


where d qx = d Tx — d rx+1 ,x = 1,... ,N — 1, and vp. z/A. We assume 
that 

MO) < K n (5.22) 

with Kn growing less than exponential in N. 

Consider a local function <p(rj) such that 

Mv)\<Cj2(pi + v (^)) S , S < 1 (5.23) 

where is the local support of ip. Let dtp be the shifted function 
(well defined for large enough N). Denoting <p(r, (3) =< <p > Tj p the 
expectation with respect to d/i^p defined by (15.81) 

Proposition 5.1. 


lim E^ 

N —^oo 


(l ds -1 


G(y)(p{T(s),p(y))dy 


= 0 . 


) 

(5.24) 


The proof uses similar ideas as in the previous cases, plus results 
contained in [13]. In particular from same argument used in [T3] follows 
that the empirical distribution for the r’s converge to the solution of 
r t r /■! i 


ds 


dyG(y, s)r(r(y, s),/3(y))dy - d y G( 1, s)r(s) 


= 0 (5.25) 


JO DO J 

for any smooth function G(y, t ) on [0,1] such that G(l) = 0 and G'(0) = 


0. 
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5.1. Some “thermodynamic” consequences. We illustrate here 
how the above limit realizes the so-called quasi-static isothermal trasfor- 
mation of the thermodynamics. When a temperature gradient is present, 
’isothermal’ should be intended that the temperature gradient does not 
change in time. In this case these are quasistatic transformations be¬ 
tween non-equilibrium stationary states. 

When performing the usual diffusive scaling, as done in [14] for con¬ 
stant temperatures and in [13] for the temperature gradient case, the 
results below are obtained in a two step limit, first the hydrodunamic 
diffusive limit, then a quasi-static limit, see details in jT3, 13] 


5.1.1. Excess Heat. The (normalized) total internal energy of the sys¬ 
tem is defined by 


N 




X=1 


pi 


+ V(r x ) 


(5.26) 


then internal energy evolves as: 

U N (t)-U N (0) = W N (t) + Q N (t) 
where ^ 

Wa r(t) = N 1+a [ f(s)p N (s)ds = [ f(s) 


dq N {s ) 


n 


is the (normalized) work done by the force r(s) up to time t, while 


N 


Q N (t) =7 N 1+a / ds (pl(s) - P x x ) 


X=1 

N 


(5.27) 


jV“ / 2 ^ v / 27^- 1 / p r ts]dv:,Uj. 


X=1 


is the total flux of energy between the system and the heat bath (di¬ 
vided by IV). As a consequence of theorem I 5. II we have that 

lim Wat (t) = ( f(s)dC(s) := W(t) 

N^oo J Q 

where C(t) = f Q l r(y, t)dy 1 the total macroscopic length at time t. While 
for the energy difference we expect that 

lim (U N (t) - U N (0)) = I [u(f(t),P(y))-u(f(0),P(y))]dy (5.28) 

n^-oo J 0 

where u(r, /3) is the average energy for f/,p tT , i.e. 

u{t,P) = I S x dp\p = jp+J V(r)e-PW-^-^dr 
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with G(t,(3) = log f e^ /3< ' V( ' r> ~ Tr) dr. Unfortunately for lack of uniform 
integrability for the energy distribution, we do not have a rigorous proof 
of (15.281) . since energy correspond to a value 5 = 2 in (15.231) . Assum¬ 
ing that the local equilibrium established in I 5.11 extends to quadratic 
growing functions, then we have obtained 

Q N (t) —► Q(t) (5.29) 

TV—>oo 

where Q(t) is deterministic and satisfy the relation 

Q(t) = [ [u{f(t),P(y))^u{f(0),P(y))}dy~yV(t). (5.30) 

Jo 

We call Q(t) excess heat (or heat in the case of constant profile of 
temperature), and (15.301) express the first principle of thermodynamics 
for quasistatic transformations. 

Notice that from (I5.27j) . Qn (£) is the time integral of a highly fluc¬ 
tuating random quantity. It is only after the particular space-time 
scaling that the this quantity converges, and that the limit is a de¬ 
terministic function of time. In this sense this result is different from 
similar identifications of heat as a stochastic flux of energy done in 
the so-called “Stochastic Thermodynamics” (cf. |T3j for example). On 
the other hand in thermodynamics heat is defined as the total flux of 
energy between the system and the thermal bath, resulting from the 
complete transformation from one stationary state to another, while it 
does not attempt to describe the instantaneous flux of energy. But in 
a quasi-static transformations at each time t a new stationary state is 
reached. 


5.1.2. Free energy and Clausius identity. Define 


F{t) = / F(r(y, t),(3(y))dx 
Jo 

A straightforward calculation gives 

J-’(t) — J-’(O) = f ds f dyf{s)d s r(y,s)= f f(s)dC(s) 
Jo Jo Jo 


(5.31) 


= m) 

(5.32) 


i.e. Clausius equality for the free energy. 

Equivalently, by using the thermodynamic relation T = u — /3~ 1 S 
from fl5.16[) . we have 

dy P~\y) ( S(r(y , t),u(y, t)) - S(r(y, 0 ),u(y, 0))) = Q(t) (5.33) 
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In the case of constant temperature profile, this reduce to the expected 
thermodynamic relation S = /3Q for quasistatic isothermal thermody¬ 
namic trasformations. 
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